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ABSTRACT 

We present a quantitative study on the properties at death of fast-rotating massive stars evolved at low- 
metallicity, objects that are proposed as likely progenitors of long-duration 7-ray bursts (LGRBs). We perform 
ID-nrotation stellar-collapse simulations on the progenitor models of Woosley & Heger (2006) and critically 
assess their potential for the formation of a black hole and a Keplerian disk (namely a collapsar) or a proto- 
magnetar We note that theoretical uncertainties in the treatment of magnetic fields and the approximate han- 
dling of rotation compromises the accuracy of stellar-evolution models. We find that only the fastest rotating 
progenitors achieve sufficient compactness for black-hole formation while the bulk of models possess a core 
density structure typical of garden-variety core-collapse supernova (SN) progenitors evolved without rotation 
and at solar metallicity. Of the models that do have sufficient compactness for black-hole formation, most of 
them also retain a large amount of angular momentum in the core, making them prone to a magneto-rotational 
explosion, therefore preferentially leaving behind a proto-magnetar. A large progenitor angular-momentum 
budget is often the sole criterion invoked in the community today to assess the suitability for producing a 
collapsar. This simplification ignores equally important considerations such as the core compactness, which 
conditions black-hole formation, the core angular momentum, which may foster a magneto-rotational explosion 
preventing black-hole formation, or the metallicity and the residual envelope mass which must be compatible 
with inferences from observed LGRB/SNe. Our study suggests that black-hole formation is non trivial, that 
there is room for accommodating both collapsars and proto-magnetars as LGRB progenitors, although proto- 
magnetars seem much more easily produced by current stellar-evolutionary models. 
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1. INTRODUCTION/CONTEXT 

All stars with masses initially between ^8 and ^150 M0 
eventually form a degenerate core that inevitably collapses 
to form a proto-neutron star. Much less certain is its subse- 
quent evolution, the potential formation of a black hole, and 
the powering of a supernova (SN) explosion, sometimes asso- 
ciated with a long-duration 7-ray burst (LGRB). The situation 
is deceptively simple and the outcome rests fundamentally on 
the solution to an energy problem. An explosion or a fizzle 
depends on the efficiency with which the system can extract 
the prodigious gravitational energy released during collapse. 
There are two principal forms of energy at disposal. The first 
one is the gravitational-binding energy liberated by the col- 
lapsing star and in particular its degenerate core. The second 
is its rotational energy (actually drawn from gravitational en- 
ergy), which is a function of the angular-momentum distri- 
bution and budget in the progenitor star. Understanding how 
these two energy sources can be channeled to power relativis- 
tic and non-relativistic ejecta in core-collapse SNe and leave 
behind a neutron star, a fast-rotating pulsar, a magnetar, or 
a black hole, has been the subject of numerous studies and 
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the source of much debate (Bethe 1990; Woosley 1993; Her- 
ant et al. 1994; Burrows et al. 1995; Janka & Muller 1996; 
Wheeler et al. 2000; Kitaura et al. 2006; Buras et al. 2006b; 
Burrows et al. 2006; Murphy & Burrows 2008; Nordhaus 
et al. 2010; Pejcha & Thompson 2011; Hanke et al. 2011; 
Takiwaki&Kotake2011). 

Thermal MeV neutrinos are abundantly radiated from the 
optically-thick, dense, and hot proto-neutron star, allowing 
its internal energy to be released on a diffusion timescale. 
In the neutrino mechanism for core-collapse SN explosions 
(Bethe & Wilson 1985), the absorption in the infalling mantle 
of ^10% of this neutrino flux may alone lead to the revival of 
the stalled shock and the ejection of the progenitor envelope 
with an asymptotic kinetic energy of up to 1 B (10^' erg) (Ki- 
taura et al. 2006; Buras et al. 2006a,b). But this generic mech- 
anism should not, however, be the origin of the larger explo- 
sion energies of ^ lOB inferred for a small fraction of core- 
collapse SNe. Instead, their scarcity calls for exceptional cir- 
cumstances, which seem intricately related to fast progenitor- 
core rotation (Burrows et al. 2007; Takiwaki & Kotake 201 1). 
It is probable that most stars contain some angular momentum 
at the time of death, either because they did not lose it com- 
pletely through the combined effects of magnetic torques and 
stellar-wind mass loss (Maeder & Meynet 2000; Meynet & 
Maeder 2000; Heger et al. 2000; Hirschi et al. 2004; Meynet 
& Maeder 2005; Hirschi et al. 2005; Heger et al. 2005), or be- 
cause they gained it from a companion star in a binary system 
(Wellstein & Langer 1999; Petrovic et al. 2005; Cantiello et al. 
2007). As the envelope collapses, the rotational energy in- 



creases. During this process, the inner core (< 0.5 Mq) spins 
up by about three orders of magnitude and remains in solid 
body rotation, while the outer core develops a differentially- 
rotating profile (Ott et al. 2006). The energy associated with 
rotation can be large (0(1 OB)) and tapped by instabilities de- 
veloping at the surface of the proto-neutron star, in particu- 
lar the magneto-rotational instability (MRI; Balbus & Haw- 
ley 1991; Akiyama et al. 2003; Obergaulinger et al. 2009). 
Numerical simulations for fast-rotating progenitor stars sug- 
gest that the magneto-rotational mechanism of explosion is 
promising and offers a very attractive explanation for the ex- 
istence of highly-energetic SNe (LeBlanc & Wilson 1970; 
Bisnovatyi-Kogan et al. 1976; Wheeler et al. 2000; Yamada 
& Sawai 2004; Moiseenko et al. 2006; Burrows et al. 2007; 
Dessart et al. 2008; Takiwaki & Kotake 2011). However, 
this mechanism relies fundamentally on the assumption that 
the MRI can increase the magnetic field to the required val- 
ues. An assumption that has not yet been shown numerically 
in the full core-collapse context, although preliminary inves- 
tigations in this sector are promising (Obergaulinger et al. 
2009). Hence, combined with the diversity of progenitor- 
core properties, these two mechanisms alone, the neutrino and 
the magneto-rotational mechanism, may explain the diversity 
of core-collapse SNe, potentially encompassing two orders of 
magnitude in explosion energy, from the low-luminosity SNe 
Il-Plateau (Pastorello et al. 2004) to highly energetic SNe Ic 
(Mazzali et al. 2002). 

A great puzzle is then to understand the necessary depar- 
tures from this general core collapse scenario to produce an 
LGRB in addition to a SN explosion, as spectroscopically 
confirmed in, to date, six LGRB/SNe pairs, (for a recent com- 
pilation, see Berger et al. 2011). This requires that ^0.1- 
1 B be injected into a low-mass, baryon-deficient collimated 
region (at the origin of 7-rays) and that at the same time 
^10 B be injected quasi-isotropically to eject the progenitor 
envelope (at the origin of the SN thermal emission observed 
in the optical). The very low occurrence rate of LGRB/SN 
per core-collapse SN of < 1% (Guetta & Delia Valle 2007) 
calls for progenitor properties that are rarely encountered in 
star formation/evolution. Interestingly, an unambiguous di- 
versity emerges among LGRB/SN observations, necessarily 
translating into a significant range for the inferred properties 
of the SN ejecta, with proposed masses and kinetic energies 
possibly varying by a factor of 5-10 for both (Berger et al. 
2011). Unfortunately, a significant uncertainty is associated 
with such inferences. For example, Iwamoto et al. (1998) 
propose an ejecta mass of 11 Mq with a total energy of 20- 
50B for GRB980425/SN1998bw, but Woosley et al. (1999) 
reproduce the light curve with an ejecta of ^5 M0 and a total 
energy of 22 B. Such differences are not surprising since both 
spectra and light curves must be modeled simultaneously and 
with allowance for the complicated non-LTE non-thermal and 
time-dependent effects controlling the radiative transfer. The 
exceptionally fast ejecta expansion of hypernovae is expected 
to strengthen the time-dependent effects seen in "standard" 
core-collapse SNe (Dessart & Hillier 2008) while the large 
production of ^*'Ni and significant mixing may sizably affect 
line-profile shapes from which the expansion rate is inferred 
(Dessart et al. 2012). 

Two LGRB central-engine models are currently favored. 
They suggest the key components for a successful LGRB/SN 
are a compact progenitor with a short light-crossing time of 
~ 1 s and fast rotation at the time of collapse. One is the col- 



lapsar model (Woosley 1993; MacFadyen & Woosley 1999; 
MacFadyen et al. 2001): A fast-rotating progenitor fails to ex- 
plode in its early post-bounce phase and instead forms a black 
hole, while the in-falling envelope eventually forms a Keple- 
rian disk feeding the hole on an accretion/viscous timescale 
comparable to that of the LGRB. It is within this disk that 
the SN explosion is triggered and the ^*'Ni synthesized. The 
other model involves a proto-magnetar (Wheeler et al. 2000; 
Bucciantini et al. 2008; Metzger 2010; Metzger et al. 2011) 
in which the LGRB is born after a successful SN explosion 
(either by the neutrino or the magneto-rotational mechanism, 
although the latter seems more likely given the rapid rotation 
required for the magnetar) and the ejection of the overlying 
envelope (or at least the onset of the ejection of the inner enve- 
lope layers that clear the proto-neutron star surface). Fast ro- 
tation in the proto-neutron star permits the huge enhancement 
of the magnetic-field energy and stresses, which strengthen as 
the proto-neutron star cools and contracts, eventually giving 
rise to relativistic ejecta (^ 10 s after the onset of collapse, 
once the neutrino driven wind decays away). In both mod- 
els, rotation is key to control the dynamics in a unique way. 
It is also key to allow the simultaneous ejection of baryon- 
deficient material at relativistic speeds over a small solid angle 
and the quasi-spherical ejection of the progenitor envelope. 

A critical difference between the collapsar model and the 
proto-magnetar model is that the collapsar has to form a black 
hole. Being so central to the model, it is legitimate to in- 
vestigate what conditions black-hole formation in this context 
(Woosley 2011). Surprisingly, little has been done on this 
problem. Numerous simulations so far have focused on the 
early collapse phase and the revival of the SN shock, stopping 
too early to make any statement concerning black-hole for- 
mation. In ID, several studies have investigated the neutrino 
signal and progenitor dependence of non-rotating failed SNe 
(Liebendorfer et al. 2004; Sumiyoshi et al. 2007; Fischer et al. 
2009; O'Connor & Ott 2011). Sekiguchi & Shibata (2011) 
have performed 2D simulations of black-hole formation and 
the subsequent formation of an accretion disk, although with 
initial conditions that are incommensurate with currently sug- 
gested LGRB progenitors. In 3D, Ott et al. (201 1) performed 
fully general-relativistic simulations of black-hole formation, 
however using a simplified soft equation of state that favored 
it. Finally, other simulations have started from a pre-existing 
black hole and investigated the powering of the jet at the ori- 
gin of the LGRB (Aloy et al. 2000; Proga et al. 2003; Zhang 
et al. 2004; Lindner et al. 2010), or the longer-term synthe- 
sis of ^''Ni in this unusual context (Milosavljevic et al. 2012). 
No simulation has ever demonstrated from first principles, and 
thus convincingly, the validity of the collapsar model, i.e., that 
the progenitors proposed for this model would indeed proceed 
through each and every necessary step: Collapse, formation 
of a proto-neutron star, failure of the shock revival, forma- 
tion of a black hole followed by that of a Keplerian disk, 
and finally the powering of both the LGRB and the SN, in- 
cluding the synthesis of a generally large amount of ^*'Ni by 
core-collapse SN standards. This is an obvious shortcoming 
of all theoretical investigations on the collapsar model and its 
proposed progenitors. Similar gaps in the modeling of proto- 
magnetar driven LGRBs exist: The early magneto-rotational 
core-collapse SN evolution has been modeled in 2D (Burrows 
et al. 2007; Takiwaki & Kotake 201 1) and so has the phase in 
which relativistic outflows are driven (Bucciantini et al. 2007, 
2008, 2009; Komissarov & Barkov 2007), but the evolution 



connecting the two phases has not been modeled. The robust- 
ness of the magneto-rotational explosion mechanism largely 
rests on the efficiency of angular-momentum transport, and in 
particular the extraction of the free-energy stored (and replen- 
ished through accretion) in differential rotation at the surface 
of the proto-neutron star The failure to extract this energy 
on short time scales may, however, facilitate black-hole for- 
mation, although it may also compromise energy extraction 
in the collapsar model. These complicated issues require de- 
tailed modeling to build upon the promising results of Thomp- 
son et al. (2005) and Obergaulinger et al. (2009). To this day, 
the collapsar model has been studied more extensively than 
the proto-magnetar model for LGRBs, so the later may look 
more promising in some ways in part because of the lesser 
scrutiny it has received. 

In this paper, we focus on one important aspect of the col- 
lapsar model to validate, or invalidate, the assumption, often 
made but so far never checked, that the LGRB progenitor 
models available in the literature indeed collapse to form a 
black hole. We do this by performing hydrodynamical simu- 
lations of the LGRB progenitor models of Woosley & Heger 
(2006) using the code GRID (O'Connor & Ott 2010). This 
issue is critical for testing the potential of progenitor stars 
for producing LGRBs via the collapsar mechanism, but may 
also serve to diagnose an attractive channel for the formation 
of proto-magnetars. Such "failed" collapsars (because they 
explode before forming a black hole) represent a serious al- 
ternative for the production of LGRBs, although they have 
their own caveats (Metzger et al. 2011). Admittedly, the phe- 
nomenon of core collapse, bounce, and the events that follow 
are fundamentally multi-D. We believe, however, that much 
can be learned from ID simulations of the kind presented 
here. For example, the mass-accretion rate onto the proto- 
neutron star is largely determined by the angular-averaged 
density profile of the progenitor star, which we capture accu- 
rately. Our ID exploration reveals the landscape of core prop- 
erties at bounce and quantifies fundamental differences be- 
tween progenitors. In the next section, we start by reviewing 
results from stellar-evolution models for LGRB progenitors. 
We then describe our methodology for the GRID simulations 
of LGRB progenitors available in the literature. We make a 
short digression in Section 4 to discuss the rotational proper- 
ties of the collapsed cores of massive stars. In Section 5, we 
present our results before concluding in Section 6. 

2. STELLAR-EVOLUTION MODELS OF LGRB PROGENITORS 

Stellar-evolution calculations have been performed to in- 
vestigate the mass, rotation, and metallicity requirements 
for producing fast-rotating pre-collapse stars. Both single- 
(Hirschi et al. 2005; Yoon & Langer 2005; Woosley & Heger 
2006; Georgy et al. 2009) and binary-star (Petrovic et al. 
2005; Cantiello et al. 2007) evolutionary scenarios have been 
investigated. Fast rotation of the proto-stellar core is clearly 
critical to procure a large angular momentum to the star ini- 
tially. If the rotation rate attained is sufficiently large, the 
star may even evolve chemically homogeneously and avoid 
a supergiant phase, which is known to sap the core of its 
rotation through the effects of magnetic torques. Such fast- 
rotating chemically-homogeneous stars also naturally die as 
H-deficient He-poor Wolf-Rayet (WR) stars. Low metallic- 
ity quenches the stellar-wind mass loss rate, a condition that 
may be more important for a single-star scenario than for 
the binary-star scenario (Yoon et al. 2006; Cantiello et al. 



2007). While the general outcome of these simulations is that 
it is possible to produce massive stars with a rapidly spinning 
core/envelope at death, it is difficult to compare the final prop- 
erties of published models. Indeed, models are rarely evolved 
all the way to an iron core. The treatments of mixing, mass 
loss, and angular-momentum loss/transport differ Magnetic- 
fields may or may not be included and when they are the pre- 
scription may differ (Spruit 2002; Zahn et al. 2007). 

Furthermore, all these studies remain speculative about the 
outcome of collapse for such progenitors. They argue for 
black hole formation and the formation of a disk based on 
order of magnitude estimates, rather than detailed numerical 
simulations. For a start, of all the LGRB progenitor simu- 
lations, only those of Woosley & Heger (2006) are evolved 
all the way to the onset of collapse. In simulations halted 
well before, the iron-core mass is estimated from the CO-core 
mass (Hirschi et al. 2005) or is simply not considered in the 
discussion (Yoon & Langer 2005). Most studies consider a 
model viable for producing a collapsar based exclusively on 
the angular- momentum budget of the inner 3M0, and per- 
form no checks on the likelihood of forming a black hole: its 
formation is deemed so obvious that the discussion of any al- 
ternate scenario is generally omitted. 

Differing in their criteria and approaches for selecting col- 
lapsar candidates, some studies may yield progenitor-mass 
ranges that do not even overlap: Using an angular-momentum 
criterion, Yoon & Langer (2005); Woosley & Heger (2006); 
Yoon et al. (2006) favor progenitor stars with a main-sequence 
mass below ^30 M©. In contrast, Hirschi et al. (2005), argu- 
ing for the need for both a large angular momentum, a large 
iron core, and a WO stellar type at death, favor progenitors 
with a main-sequence masses above ^35 M0 (magnetic fields 
are not treated in this work, though). 

Recent studies suggest that selecting collapsar progenitors 
based exclusively on a large angular-momentum budget may 
be too simplistic. For example, the magneto-rotational explo- 
sion invoked to explain hypernovae derives its energy from 
this same large core angular momentum (via the MRI and 
the strongly differentially-rotating layers in the post-shock re- 
gion). This mechanism does not obviously accommodate the 
formation of a black hole, as demonstrated by Dessart et al. 
(2008). They simulated the collapse of the core and the de- 
velopment of a magneto-rotational explosion in model 350C 
of Woosley & Heger (2006) and found that rotational energy 
of order 10 B is readily available to launch a SN ejecta on 
a timescale of a few 100 ms. They furthermore found that 
accretion is easily shut off by the developing explosion and 
that the proto-neutron star mass fails to grow to the instability 
threshold for black-hole formation. However, the simulations 
of Dessart et al. (2008) did not resolve the MRI but instead 
used an equipartition ansatz to estimate the magnitude of the 
MRI-amplified magnetic fields. In reality, magnetic-field re- 
connection may, for example, compromise the dynamical po- 
tential of magnetic stresses, channeling magnetic energy into 
heat to be radiated away by neutrinos. This and other alterna- 
tives have been studied by Thompson et al. (2005) under the 
general form of viscous dissipation. They find that the extra 
energy deposition can, in some cases, considerably alter the 
post-bounce dynamics and generate a vigorous explosion. 

O'Connor & Ott (201 1) studied black-hole formation based 
on a variety of progenitor models characterized by differ- 
ent main-sequence mass, metallicity, and rotation rate. They 



find that the outcome of collapse can be anticipated from the 
bounce compactness of the progenitor core structure and in 
particular that of the region inside 2.5 Mq, which corresponds 
approximately to the maximum mass that a proto-neutron star 
can have and remain in hydrostatic equilibrium. They find 
that higher mass progenitors published in the literature do 
not always have larger iron cores and therefore that they are 
not necessarily more prone to black-hole formation. They 
also reveal considerable diversity in progenitor core structure, 
even for the same main-sequence mass. Some stellar evo- 
lution studies obtain a monotonic increase of the iron-core 
mass (or bounce compactness; see Fig. 9 of O'Connor & Ott 
2011) versus main-sequence mass (Limongi & Chieffi 2006; 
see also Hirschi et al. 2004, 2005), while some show an anti- 
correlation beyond ^40 M0 (Woosley & Heger 2007). The 
primary reason for this are differing prescriptions for rate and 
time of mass loss, one of the major uncertainties in massive 
star evolution (see also the discussions in Hirschi et al. 2005 
and O'Connor & Ott 2011). 

3. METHODS & INITIAL MODEL SET 

In this work, we use the open-source, spherically symmet- 
ric, general relativistic, Eulerian hydrodynamics code GRID 
(O'Connor & Ott 2010). Rotation is included through a 
centrifugal-acceleration term in the momentum equation — 
this is the most important dynamical feature of rotation rel- 
evant to core collapse. However, GRID cannot account for 
the associated deviations from spherical symmetry nor any 
angular-momentum redistribution. We select the equation of 
state (EOS) from Lattimer & Swesty (1991) characterized by 
a nuclear incompressibility of 220 Me V (hereafter referred to 
as the LS220 EOS). This EOS provides the best match to both 
mass and mass-radius constraints from observations and nu- 
clear theory (Demorest et al. 2010; Ozel et al. 2010; Steiner 
et al. 2010; Hebeler et al. 2010). GRID uses an efficient neu- 
trino leakage/heating scheme that qualitatively reproduces the 
salient features of neutrino transport. We refer the reader to 
O'Connor & Ott (2010, 201 1) for additional details on GRID 
and our methodology. 

As described above, the only stellar-evolutionary models 
for LGRB progenitors that are evolved until the onset of col- 
lapse are those proposed by Woosley & Heger (2006). We 
thus focus on their model dataset for our investigation on the 
dynamics of the core-collapse SN engine and the potential 
formation of a black hole in the collapsar context. Using 
KEPLER, Woosley & Heger (2006) investigated a rather nar- 
row range of progenitor masses, but varied the initial rotation 
rate (solid-body rotation is assumed initially) and environ- 
mental metallicity from solar to 1% solar (with an additional 
tunable factor as low as 0. 1 for the metallicity-dependent mass 
loss rate, equivalent to a reduction in metallicity by a factor of 
100 in their mass-loss prescription). Arguing that the inferred 
mass of LGRB/SN ejecta known in 2006 is on the order of 
IOM0, and since higher-mass stars may lose too much angu- 
lar momentum through stellar winds (even at low metallicity), 
they focused primarily on lower-mass progenitors, with main- 
sequence masses of 12 and 16 M©,^ with the exception of one 
35 Mq model set. 

We adopt the same nomenclature as for their 12-, 16- and 
35-Mq models. It comprises the model's main-sequence 

' They also perform simulations for I6-M0 helium cores and find compa- 
rable outcomes. 
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Figure 1. Angular velocity f2(r) versus radius r at both the pre-SN stage 
(dashed lines) and at core bounce (solid lines) for selected models of Woosley 
& Heger (2006). The inner homologously collapsing core maintains its initial 
uniform rotation throughout collapse. 

mass, followed by a letter denoting the environmental metal- 
Hcity ('S' for solai-, 'O' for 10% solar, and 'T' for 1% solar). 
An additional letter is appended to individualize the models 
done with different WR mass-loss rate prescriptions, allow- 
ing or not for magnetic effects, and the total angular momen- 
tum of the star I6-M0 helium models are denoted by 'HE16' 
followed by an individualizing capital letter. 

In this work, we simulate the collapse and post-bounce evo- 
lution with GRID for all these progenitor models, with a pri- 
mary focus on determining their ability to produce the key 
features of the collapsar model: A black hole together with a 
Keplerian disk. As we discuss in the following section, black- 
hole formation is not obviously guaranteed in any of these 
dying stars. 

4. NOTES ON ROTATING CORE COLLAPSE 

Since LGRBs seem fundamentally related to rapid rotation, 
it is useful to summarize a few facts and concepts related to 
the gravitational collapse of rotating iron cores in massive 
stars. First, it is reasonable to assume (which is borne out by 
simulations, e.g., Heger et al. 2005) that the iron core, in its 
pre-collapse state, will be approximately uniformly rotating. 
Such a solid-body rotation corresponds to the lowest energy 
state at fixed total angular momentum and will be assumed on 
a secular timescale by any rotating fluid that has some means 
to redistribute angular momentum. 

Rotating core collapse, even for the high pre-collapse rota- 
tion rates of some of the potential LGRB progenitors that we 
consider in this study, proceeds qualitatively in a very simi- 
lar fashion to non-rotating collapse as long as the ratio of the 
centrifugal acceleration Ocent to the gravitational acceleration 
flci-av, is small. 
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Due to angular-momentum conservation, the angular velocity 
behaves as ^{f) ex r"^. M{r) stays constant for a collapsing 
mass shell, and, thus, the above ratio increases during collapse 
as r"' and may potentially become large for small radii. 
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0.182 


24.9 


3.43 


0.007 


12TH 


0.107 


1.042 






(1.93) 


(1.79) 


6.67 


0.495 


85.2 


2.56 


0.138 


5.06 


3.59 


0.313 


12TI 


0.145 


1.323 






(2.02) 


(1.86) 


3.33 


0.507 


15.3 


7.46 


0.144 


3.26 


3.48 


0.610 


12TJ 


0.517 


1.281 


0.853 


2.51 


2.37 


2.97 


0.640 


2.64 


8.58 


0.191 


1.10 


4.31 


3.286 


16SA 


0.101 


0.000 




(1.88) 


(1.74) 


(14.6) 








0.138 








16SG 


0.109 


0.203 




(1.91) 


(1.77) 


9.79 


0.404 


2.71x10'' 


2.16 


0.141 


20.8 


3.27 


0.010 


16SH 


0.081 


0.341 




(1.76) 


(1.64) 


(7.70) 








0.182 


16.8 


2.68 


0.019 


16SI 


0.380 


0.189 


1.132 


2.38 


2.20 


(9.85) 








0.158 


11.9 


3.82 


0.062 


16SL 


0.075 


0.207 




(1.73) 


(1.62) 


(6.30) 








0.130 


28.4 


3.31 


0.007 


16SM 


0.121 


0.229 




(2.02) 


(1.85) 


(8.31) 








0.145 


20.0 


3.51 


0.016 


16SN 


0.496 


0.455 


0.777 


2.42 


2.27 


9.45 


0.508 


40.6 


1.77 


0.187 


3.25 


3.69 


0.451 


160A 


0.144 


0.000 




(2.15) 


(1.96) 


(15.8) 








0.133 








160G 


0.193 


0.176 


3.437 


2.33 


2.11 


7.16 


0.230 


4.37x10"^ 


8.49 


0.168 


17.6 


3.82 


0.018 


160H 


0.185 


0.248 




(2.21) 


(2.00) 


9.18 


0.133 


810. 


7.9x10-5 


0.150 


20.5 


3.59 


0.023 


1601 


0.344 


0.733 


1.449 


2.37 


2.19 


7.10 


0.553 


18.2 


5.11 


0.152 


3.21 


4.06 


0.700 


160L 


0.124 


0.316 




(2.02) 


(1.86) 


8.63 


0.200 


593. 


0.05 


0.138 


14.1 


3.41 


0.030 


160M 


0.172 


1.059 




(2.17) 


(1.98) 


5.64 


0.590 


24.5 


6.31 


0.177 


2.60 


3.02 


0.480 


160N 


0.357 


1.382 


1.458 


2.40 


2.22 


3.36 


0.582 


4.70 


10.8 


0.162 


1.40 


3.59 


2.082 


16TA 


0.070 


0.000 




(1.76) 


(1.64) 


(16.0) 








0.148 








16TG 


0.288 


0.242 


1.738 


2.35 


2.16 


13.3 


0.366 


3.44X 10"* 


2.41 


0.174 


10.4 


4.02 


0.083 


16TH 


0.434 


0.598 


0.958 


2.41 


2.25 


8.01 


0.511 


23.7 


3.44 


0.151 


2.57 


3.80 


0.599 


16TI 


0.242 


1.367 


2.791 


2.41 


2.21 


3.51 


0.554 


10.6 


10.4 


0.150 


2.17 


3.77 


1.341 


350A 


0.178 


0.289 




(2.26) 


(2.05) 


(12.9) 








0.153 


14.6 


3.53 


0.044 


350B 


0.537 


1.545 


0.776 


2.42 


2.25 


16.4 


0.545 


31.5 


4.80 


0.198 


1.44 


3.48 


3.617 


350C 


0.458 


1.980 


0.972 


2.43 


2.29 


4.44 


0.622 


4.84 


23.6 


0.162 


1.10 


4.49 


7.521 


HE16C 


0.137 


0.133 




(2.06) 


(1.90) 


(5.15) 








0.134 


25.5 


3.11 


0.007 


HE16D 


0.283 


0.440 


1.706 


2.35 


2.16 


8.65 


0.367 


117. 


0.88 


0.171 


4.74 


3.65 


0.206 


HE16E 


0.129 


1.428 




(2.00) 


(1.85) 


6.82 


0.594 


43.1 


6.05 


0.153 


3.77 


3.27 


0.447 


HE16F 


0.496 


1.096 


0.837 


2.48 


2.33 


4.29 


0.567 


8.97 


10.5 


0.202 


1.46 


4.28 


2.323 


HE16H 


0.610 


1.196 


0.641 


2.52 


2.38 


4.12 


0.597 


5.95 


11.6 


0.204 


1.26 


4.23 


3.335 


HE16K 


0.132 


0.134 




(2.04) 


(1.88) 


(5.16) 








0.140 


26.6 


3.11 


0.007 


HE16L 


0.316 


0.315 


1.497 


2.36 


2.17 


9.34 


0.286 


195. 


0.24 


0.165 


6.73 


3.73 


0.116 


HE16M 


0.111 


1.206 




(1.91) 


(1.77) 


10.4 


0.532 


94.0 


2.61 


0.134 


4.42 


3.23 


0.348 


HE16N 


0.198 


1.203 


3.424 


2.36 


2.15 


7.81 


0.604 


35.9 


7.15 


0.154 


2.59 


3.84 


0.970 


HE160 


0.298 


1.209 


1.891 


2.39 


2.20 


6.77 


0.594 


25.0 


8.86 


0.127 


1.98 


3.61 


1.414 


HE16P 


0.573 


1.038 


0.( 


)72 


2.49 


2.35 


6.42 


0.584 


14.2 


9.46 


0.235 


1.58 


4.16 


2.523 



^ Model designation from Woosley & Heger (2006). See text for details. 

Initial central angular velocity. 
'^ Time elapsed between bounce and black hole formation, ■ ■ ■ indicates that no black hole formed within 3.5 s of bounce. 

Baryonic mass of the proto-neutron star at the time of black hole formation. If no black hole forms in 3.5 s, we give the proto-neutron star baryonic mass at 3.5 s. 
^ Gravitational mass of the proto-neutron star at the time of black hole formation. If no black hole forms in 3.5 s, we give the proto-neutron star gravitational mass at 3.5 s. 

Baryonic mass interior to the innermost stable circular orbitat the time of disk formation. If the angular momentum is too low to foster disk formation, we give the progenitor mass 
in parentheses instead. 
^ Dimensionless spin of the black hole when the disk forms. ■ ■ ■ indicates no disk forms. 

Twice the free-fall time of the mass element at the innermost stable circular orbit. ■ ■ ■ indicates no disk forms. 
^ Baryonic mass outside of the black hole at disk formation. ■ ■ ■ indicates no disk forms. 

J Rotational period Prcf = 2-k ju) computed by appoximating u) as the ratio of the total angular momentum to the moment of inertia of the proto-neutron star, Eq. 1 . ■ ■ ■ indicates a 
non-rotating model. 
" Proto-neutron star moment of inertia, note this can be up to two times the value for a non-rotating cold neutron star. This is the origin of the discrepancy with the values presented 



by Woosley & Heger (2006), who consider a cold non-rotating 1.4Mq neutron star with a moment of inertia/ = 1.4x 10 



gem . 



indicates a non-rotating model. 



Free energy stored in differential rotation. This amounts to the energy difference between the proto-neutron star we obtain with GRID and the corresponding proto-neutron star with 
the same total angular momentum and moment of inertia but assuming solid-body rotation. ■ ■ ■ indicates a non-rotating model. 



In the case of Ocent/flgrav ^ 1, the collapsing rotating iron 
core will behave like a non-rotating core and separate into a 
subsonically collapsing inner core (|vr(r)| < Cs(r)) and a su- 
personically collapsing outer core (|vr(r)| > Csir)). The in- 
ner core exhibits a self-similar (homologous) velocity profile, 
v(r) ex r, until core bounce and shock formation (Goldreich 
& Weber 1980). After core bounce, the inner core material 
forms the core of the proto-neutron star and outer core ma- 
terial accumulates at its edge. The mass of the inner core at 
bounce is typically ^Q.SMq for non-rotating cores, set by the 
EOS of relativistic electrons, the trapped lepton number, core 
entropy, and gravity (Burrows & Lattimer 1983). It increases 
monotonically, though slowly, with increasing pre-collapse 
rotation rate. For cores that reach flcent/ograv ~ 1, the mass of 
the inner core will be increased to > 0.7 Mq (Dimmelmeier 
et al. 2008). 

Since the inner core is collapsing homologously, we can 
introduce a homology parameter a(t), so that 



r(t) = a(t)rQ 



(2) 



where r{t) is the radius of a collapsing fluid element at time t 
and ro is its initial radius. This must hold for any mass shell 
within the inner core. Two collapsing fluid elements located 
initially at ri and r2, and rotating with a frequency ili = j\/r\ 
and Vt2 = ji/r^ conserve angular momentum, Q,\ rj = ][ = j\ = 
r^iTp Homology implies r[ = a{t)ri, and r!^ = a{t)r2, therefore 



^('■'i) _ ii/'-'i _ ii ««' 



^li'-'i) jilr'2 h o^it) 



2„2 






(3) 



Since this property holds for any mass shell within the inner 
core, the rotational profile must be preserved under homol- 
ogous collapse. In practice, because the pre-collapse cores 
of massive stars are always in solid-body rotation, so are the 
inner cores of the proto-neutron star at bounce. Outside of 
the homologously collapsing core the self similar relation of 
Eq. 2 does not hold. Most generally, gradients in the rota- 
tion rate develop due to the underlying density gradients in 
the hot postshock region and in the supersonically infalling 
region ahead of the accretion shock. 

Based on these arguments, we expect an early post-bounce 
rotational profile that is approximately uniform within the in- 
ner O.5-O.7M0 (out to ^10- 15 km in radius) and strongly 
differential at larger radii. This is confirmed by Fig. 1 which 
shows n(r) at bounce as obtained with GRID for a variety of 
models considered in this study. We also show the rotational 
profile at the onset of collapse. This result is not entirely new, 
but has previously been pointed out by Ott et al. (2006) in the 
context of ID and 2D rotating core collapse simulations. 

Since uniform rotation is the lowest energy state, the shear 
energy of differential rotation is to be interpreted as a free en- 
ergy that will be tapped by any process (e. g., nonaxisymmet- 
ric rotational shear instabilities, viscosity, or the MRI) capable 
of redistributing angular momentum. Viscosity would lead 
to additional heating in the postshock region to enhance the 
neutrino mechanism (Thompson et al. 2005) while the MRI 
action could strengthen the magnetic fields, driving bipolar 
outflows in the magneto-rotational mechanism (Burrows et al. 
2007). 

In our simulations, we estimate the available free energy of 
differential rotation by computing the difference in rotational 



energy of the proto-neutron star model in GRID and the ro- 
tational energy of a uniformly spinning proto-neutron star of 
the same angular momentum and moment of inertia. 



Frn, = T- 






phXW^lr^dr , 



where from O'Connor & Ott (2010), for GRID, 
87r 

rRpNs I pRpNs 



1= 



^PNS 



phXW^r^dr , 



(4) 



(5) 



(6) 



UJ = 



phXW^rv^r^dr / / phXW^r'^dr (7) 







and T is the rotational energy, / is the moment of inertia, and 
cj is the uniform rotation frequency, h is the specific enthalpy, 
X^ is the grr component of the metric, W is the Lorentz factor, 
and Vip is the angular velocity. We take /?pns to be the radius 
where the matter density, p = 10'*'gcm"-'. 

5. RESULTS 

We have performed simulations for the entire set of 
KEPLER models published in Woosley & Heger (2006)*". We 
first consider the rapidly spinning progenitors evolved with- 
out magnetic fields. All these models have a dimensionless 
Kerr spin (a* = Jc/GM^) at 3 M© greater than unity (with 
the exception of model HE16J, which has a* = 0.91) and are 
thus considered as promising collapsar candidates by Woosley 
& Heger (2006). Unfortunately, when evolved with GRID, 
the collapsing iron core of all such models halts its collapse 
and expands — these models do not experience core bounce 
within a few seconds of evolution in GRID. We associate this 
problem with the neglect of the centrifugal acceleration in the 
momentum equation in KEPLER, an approximation that fails 
in the fastest rotating models. This term is included in GRID. 
The mismatch suggests that their fastest models may be sig- 
nificantly affected by the addition of this term. Even if they 
did collapse, it is not clear that such extremely fast rotating 
cores would avoid a centrifugal bounce. 

In the remainder of this paper, we thus limit our discussion 
to models evolved with magnetic fields and therefore subject 
to magnetic torques during their evolution. Of the 46 mod- 
els that fulfill this criterion, we identify 4 additional models 
(120M, 16TJ, 350D, and HE16G) that do not collapse but in- 
stead expand when restarted with GRID. We exclude these as 
well from our study. Finally, for reference and completeness, 
we include the non-rotating models associated with each se- 
ries (12SA, 120A, 12TA, 16SA, 160A, 16TA), making a total 
of 48 models. Each simulation is continued after core bounce 
until a black hole forms or until a time of 3.5 s has passed, 
whichever comes first. We present the results for these 48 
models evolved with magnetic fields in Table 1 (for the table 
layout, we group models in bundles first of increasing mass, 
then of decreasing metallicity, and finally in alphabetical or- 
der which generally corresponds to an increased initial rota- 
tion rate). 

As advocated by O'Connor & Ott (2011), the bounce com- 
pactness is a robust quantity for diagnosing the propensity to 

* Models are available from http://homepages.spa.umn.edu/~alex/GRB2/ 



black-hole formation, which is largely determined by the spa- 
tial extent encompassed by the 2.5 M© Lagrangian mass co- 
ordinate at core bounce in the progenitor core. To avoid in- 
troducing biases associated with the non-uniform conditions 
in the progenitor simulations (KEPLER models are not all 
evolved to the same central density on their collapse trajec- 
tory) this compactness is unambiguously evaluated at the time 
of bounce. The formal definition of this core compactness is 



?M- 



M/Mc^ 



^(Mtery=M)/1000km 



'='bo, 



(8) 



where we take M = 2.5Mq. ^(Mbary = 2.5Mq) is the radial 
coordinate that encloses 2.5 M© of baryonic material at the 
time of core bounce (O'Connor & Ott 201 1). 

Our simulations first demonstrate that most of the mod- 
els have a small core compactness ^2.5- O'Connor & Ott 
(201 1) argues that a compactness of 0.45 represents a thresh- 
old value, for the neutrino mechanism, since above it an unre- 
alistic neutrino-heating efficiency is required to prevent black- 
hole formation. We further confirm this by determining the 
critical heating efficiency for the models in Table 1 via the 
same procedure as in O'Connor & Ott (2011), to which we 
refer the reader for full details. We note that this criterion 
is for explosions via the neutrino mechanism and therefore 
neglects any magneto-rotational contribution to the powering 
of an explosion (Dessart et al. 2008), so this threshold value 
is probably a lower limit. We also stress that this criterion 
is based on spherically symmetric simulations with an effi- 
cient, but crude, approximation to neutrino transport. There- 
fore we suggest caution when attempting to interpret the out- 
come of a model based solely on the bounce compactness, 
since the adopted threshold value for BH formation is a semi- 
quantitative estimate based on where the slope of the required 
heating efficiency begins to clearly increase (see the discus- 
sion in O'Connor & Ott 2011). When plotting the critical 
heating efficiency of the non-rotating solar metallicity stars 
from Woosley & Heger (2007) determined in O'Connor & 
Ott (2011), together with the generally fast-rotating progen- 
itors of Woosley & Heger (2006), we find that both datasets 
in fact overlap for the most part (Fig. 2). In other words, in 
terms of compactness, most of these progenitors are similar 
to garden variety, low-mass, non-rotating, progenitors and do 
not seem to have any more reason to form a black hole than, 
e.g., the RSG progenitors expected to produce SNe Il-Plateau. 
As shown in Table 1 , provided no explosion is launched, half 
of these models have not formed a black hole after 3.5 s, and 
only ^15% do within ~'l s. We note that this result is not so 
surprising given the small iron-core mass (^ \AMq) of most 
Woosley & Heger (2006) models (see their Tables 1 and 2). 
The first conclusion from our work is therefore that most of 
the models presented here are rather unlikely to form a black 
hole and thus may fail in a very fundamental way to produce 
a coUapsar, urespective of their angular-momentum budget. 

Within each sequence presented in Table 1 , the models that 
form a black hole within 3.5 s of core bounce, and thus at least 
in principle susceptible to form a coUapsar, are the faster ro- 
tating ones characterized by very weak mass-loss rates. These 
properties conspire to produce larger CO cores, more typical 
of more massive stars that do not rotate. In the following dis- 
cussion, we group these models into several categories. 

The first category are models which obviously do not give 
rise to a LGRB, either by the standard Type-I collapsar or 
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Figure 2. Illustration of the critical heating efficiency ti"^1^^ versus bounce 
compactness §2.5 for our GRID simulations of the Woosley & Heger (2006) 
models, whose properties are summarized in Table 1 (blue diamonds). For 
comparison, we overplot the same quantity for the standard non-rotating core- 
collapse SN progenitor models of Woosley & Heger (2007) evolved at solar 
metallicity (red squares). For the most part, the two distributions overlap, 
suggesting that the propensity to black-hole formation and explosion is com- 
parable for both. Only models with the fastest rotation rates achieve a larger 
compactness in excess of 0.4-0.5, but these may then be diverted from black- 
hole formation through an early magneto-rotational explosion. 

the proto-magnetar mechanism, because they contain too little 
angular momentum. Optimistically assuming a failed core- 
collapse SN, which is unlikely given the modest values of 
^2 5, models 12SG (6 5 = 0.239), 160G (^s = 0.193), 16SI 
(6.5 = 0.380), and 16TG (C2.5 = 0.288) possess too little an- 
gular momentum in the remainder of the star to form a disk 
about the central black hole within 10^ s of collapse. This 
behavior is reflected by the stellar type at the time of death, 
i.e. a BSG star for model 12SG and a RSG star for models 
160G and 16TG, only 16SI is a WR star at the time of death. 
Quantitatively, this can be further inspected in Table 1 where 
we include the disk formation time, the black hole mass and 
spin at that time, and the mass exterior to the disk. We define 
disk formation to be when the accreting material will first be 
supported at the innermost stable circular orbit about a black 
hole with the enclosed mass and angular momentum using the 
formulae of Bardeen et al. (1972). We estimate the disk for- 
mation time as twice the free fall time of the innermost mass 
element that reaches a Keplerian velocity (O'Connor & Ott 
2011;BuiTows 1986). 



fDF = 2x t:\ 



'[Ve-SN(Mdi.sk)]^ 

8GMai.w 



(9) 



where rp,e-sN(A^di,sk) is the radius of the disk-forming La- 
grangian mass element in the pre-SN model. If no such mass 
element exists, no disk will form. In this case we include, 
instead of the enclosed black-hole mass, the total pre-SN stel- 
lar mass in parentheses. In the four models mentioned above, 
either no disk forms or the disk formation time is > 10^ s. 

Additionally, we can discuss the potential for these models 
to form a LGRB via the proto-magnetar model. Using Eqs. 4 - 
7, we calculate the free energy available in differential rota- 
tion at 100 ms after bounce. We also calculate a reference spin 
period (Pref = Ztt/w), measured at the onset of the neutrino 
driven explosion, by assuming solid-body rotation for the en- 
tire proto-neutron star with the same total angular momentum 



8 



and moment of inertia. The corresponding values are given in 
Table 1. In Fig. 3, we show the free energy available in dif- 
ferential rotation at 100 ms and the spin period of the proto- 
neutron star at the onset of explosion. The total rotational 
energy of the proto-neutron star, estimated as luP' /2, will in- 
crease as the proto-neutron star cools and contracts. In mod- 
els 12SG, 160G, 16SI, and 16TG, which are contained within 
the green (lightest shade) box of Fig. 3, <0.1B of free energy 
could be extracted from differential rotation via the MRI and 
converted to explosion energy, much less than is needed for 
a magneto-rotational explosion. Also, the proto-neutron star 
spin periods are > 10 ms,^ significantly larger than the < 2 ms 
periods required for the proto-magnetar model to reproduce 
classical LGRB energies (Metzger et al. 2011). 
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Figure 3. Reference proto-neutron star spin period, P|-ef, taken at the onset 
of explosion (left axis, blue dots; Eq. 7) and the free energy stored in differ- 
ential rotation 100 ms after bounce fjj,?"""' (right axis, red stars; Eq. 4) versus 
bounce compactness ^2.5 f™ all rotating models in Table 1. While models 
with a low bounce compactness show a diversity in core-rotation proper- 
ties, those with a high bounce compactness systematically have short spin 
periods and a large budget of free energy stored in the differential rotation. 
Shaded boxes refer to specific groupings of models discussed in the text. Us- 
ing §2.5 > 0.45 as a black-hole formation criterion for non-rotating progen- 
itors, we can qualitatively compare the reference spin periods of this figure 
to Metzger et al. (2011), who sketches the outcome of collapse as a function 
of progenitor spin and mass. From this, one would predict that none of the 
LGRB progenitor models studied here formed black holes. 

The second category are models with a larger angular- 
momentum budget but unfavorable bounce compactness. Al- 
though compact enough to lead to black-hole formation 
within 3.5 s of core bounce, we find that the predicted crit- 
ical heating efficiencies are similar to that expected for a 
standard 15 M© non-rotating RSG progenitor star (O'Connor 
& Ott 2011). These properties make them unlikely collap- 
sar progenitors, but in contrast, make them ideal candidates 
for proto-magnetar formation, and perhaps LGRBs through 
that channel. These models include 160I (^25 = 0.344), 
160N (a 5 = 0.357), 16TH (^ 1 = 0.434), 16TI (a 5 = 0.242), 
HE16D (a 5 =0.283), HE16L (^ 5 =0.316), HE16N (^5 = 
0.198), and HE160 (6.5 = 0.298) and are contained in the or- 
ange (medium shade) box of Fig. 3. In addition to having crit- 
ical heating efficiencies similar to what is needed to explode 

' Even taking into account the spin up due to the PNS cooling and contrac- 
tion, which will decrease the moment of inertia from the value in Table 1 to 
~ 0.4MpNsRpNs ~ 1-6 X 1O'*^(M/1.4M0)(R/I2km)- (Metzger et al. 201 1), 
or roughly a factor of 2, the spin periods are > 5 ms. 



typical low-mass massive stars, the free energy available in 
rotation is 0(1 B). This energy may be converted to explosion 
energy via the magneto-rotational mechanism. The spin pe- 
riod of these proto-neutron stars is in the range 1-6 ms, thus 
on the order of what is needed for the proto-magnetar model 
of LGRBs (Metzger et al. 201 1). 

Eventually, the fastest rotating progenitor models evolved 
with a strongly inhibited stellar-wind mass loss represent 
more suitable collapsar candidates, although each model has 
caveats. This set is contained in the purple (darkest shade) 
box of Fig. 3 and includes models 12TJ (^5 = 0.517), 16SN 
(a.5 = 0.496), 350B (^.5 = 0.537), and 350C (6.5 = 0.458). 
Model 12TJ will form a 2.37 M© (gravitational mass) black 
hole 0.85 s after core bounce, followed by a Keplerian disk 
after 2.64s, with a potential ejecta mass of 8.57 Mq. How- 
ever, much like the models in the previous category, model 
12TJ has ^3 B of free energy available in rotation that may 
lead to a magneto-rotational explosion early-on, preventing 
collapsar formation. This model is evolved at 1 % solar metal- 
licity, with an additional mass-loss rate scaling of 0.1, equiv- 
alent to an overall evolution at 10"^ solar metallicity, much 
below that observed for LGRB/SN sites. We find that mod- 
els HE16F, HE16H, and HE16P have similar characteristics 
to model 12TJ. Model 16SN forms a 2.27 M© black hole 
0.78 s after bounce. Being evolved at an effective metallic- 
ity of 0.01 solar, it has a lower angular-momentum budget at 
death and is thus more likely to avoid a magneto-rotational 
explosion. However, it forms a Keplerian disk only 40.6 s 
after core bounce, with only 1.77 M© left over for the SN 
ejecta. Such characteristics might in fact be more amenable 
to reproduce recent observations of LGRB/SNe characterized 
by a very early and narrow light-curve peak, as witnessed 
for example for GRB100316D/SN 2010bh (Chornock et al. 
2010). They may even explain why no SN is found in as- 
sociation with some nearby LGRBs (Fynbo et al. 2006). Fi- 
nally, models 350B and 350C form a black hole within 0.78 
and 0.97 s of bounce, respectively. Model 350B will accrete 
^ 16.4 M© before a Keplerian disk forms ^31.5 s after the 
onset of collapse, 4.8 M© is then available for the SN ejecta. 
With the 350C model, a disk forms very quickly after col- 
lapse, in 4.8 s, and a significant amount of mass is exterior to 
the disk, 23. 6M©, and thus much too large to accommodate 
inferred LGRB/SN ejecta masses. However, the propensity 
to collapsar formation of the 350B and 350C model may be 
ill-founded if the MRI is successful at powering a magneto- 
rotational explosion. The free energy available in rotation is 
huge, i.e. on the order of 4-7.5 B. In fact, in the 2D magneto- 
hydrodynamic simulations of Dessart et al. (2008) based on 
the 350C model, it was found that, despite the large progeni- 
tor compactness, a magneto-rotational explosion was initiated 
^200 ms after core bounce and that the proto-neutron star 
mass decreased thereafter, never reaching the mass thresh- 
old for black-hole formation. In our models, the protoneu- 
ti-on stars in models 350B and 350C have -- 30-70 B of to- 
tal rotational energy at the onset of explosion, amply match- 
ing the inferred energies of observed hypernovae. These in- 
ferences are based on the assumption that energy extraction 
from the differentially-rotating layers at the proto-neutron star 
surface is efficient and can power an explosion. Failing to 
do so, black-hole formation would result, although the ques- 
tion of energy extraction from the disk for the powering of 
a GRB would then arise. Detailed multi-dimensional core- 
collapse simulations need to be carried out to investigate the 



9 



efficiencies of magnetic/rotational/hydrodynamical instabili- 
ties for the transport of angular momentum and the extraction 
of rotational energy. 

6. DISCUSSION 

In this paper, we have performed ID general-relativistic hy- 
drodynamical simulations with GRID of the collapse, bounce, 
and post-bounce phases of the LGRB candidates of Woosley 
& Heger (2006) to investigate their propensity to black-hole 
and disk formation. We find that these progenitors are at odds 
with the proposed criteria for a collapsar progenitor or with 
the inferred properties of observed LGRB/SNe, namely a H- 
deficient He-poor WR star with a massive iron core (equiva- 
lent to a large compactness), a large angular momentum to 
form a Keplerian disk soon after black-hole formation, an 
ejecta of ~10Mq, and an evolution at about 0.1 solar metal- 
licity. A critical aspect that we focus on in this study is the 
compactness of the progenitor cores at bounce, a quantity that 
helps diagnose the likelihood of black-hole formation. 

We group the Woosley & Heger (2006) models in different 
categories according to their suitability for producing collap- 



1. Models with a dimensionless Kerr spin parameter 
greater than unity at an enclosed mass of 3M0, i.e. 
the models identified by Woosley & Heger (2006) as 
having the best potential for collapsar formation, fail 
to collapse when evolved with GRID. Their cores are 
so fast spinning that the associated centrifugal acceler- 
ation leads them into expansion. All these models are 
evolved until death in KEPLER without magnetic fields 
and centrifugal forces, which seems questionable given 
the unrealistically short spin periods at collapse. 

2. Models evolved with magnetic fields produce much 
lower rotation rates and most collapse with GRID. 

3. Of those evolved with magnetic fields, models with 
moderate rotation produce progenitors with a small 
compactness comparable to that characterizing the low- 
mass massive-star models proposed as progenitors of 
garden-variety core-collapse SNe. A small fraction of 
these is endowed with sufficient angular momentum to 
make a proto-magnetar, and thus a potential channel for 
producing LGRBs. 

4. A few models (12TJ, 16SN) with the fastest rotation 
possess a large compactness favorable for black-hole 
formation and sufficient angular momentum for the for- 
mation of a Keplerian disk, but they require evolution at 
metallicities in the range 0.0001-0.01 Zsoi, significantly 
lower than the metallicity of a few tenths solar or even 
higher at which these LGRB/SNe are found (Modjaz 
et al. 2008; Levesque et al. 2010b, a). Although in many 
respects very attractive for forming a collapsar (if we 
ignore its huge core angular momentum), model 350C 
is characterized by a large envelope mass of ^23 M0, 
which is a factor 2-10 times larger than the inferred 
ejecta mass of LGRB/SNe discovered so far (for a sum- 
mary, see Berger et al. 2011). We note that these models 
have a large angular momentum in the core, as mod- 
els in the previous category, and may thus experience a 
magneto-rotational explosion preventing collapsar for- 
mation. 



Our quantitative study thus spells out the various shortcom- 
ings of these progenitor stars for producing collapsars. Even 
in those models that have the right compactness for black- 
hole formation and sufficient angular momentum for disk for- 
mation, it is still unresolved today how they would avoid 
the magneto-rotational mechanism of explosion that is used 
to explain hypernovae (LeBlanc & Wilson 1970; Bisnovatyi- 
Kogan et al. 1976; Wheeler et al. 2000; Yamada & Sawai 
2004; Moiseenko et al. 2006; Burrows et al. 2007; Dessart 
et al. 2008; Takiwaki & Kotake 201 1) The difficulty of form- 
ing a black hole and avoiding a magneto-rotational explosion 
in fast-rotating cores, at least in the models of Woosley & 
Heger (2006), lends credence to the proto-magnetar model of 
LGRB/SNe. 

Uncertainties in mass loss at low metallicity, and in partic- 
ular during transient phases of dynamical mass loss as ob- 
served in some Luminous Blue Variable stars, is an issue, 
since it may completely dominate the mass lost in the form of 
a weaker, but secular, steady-state wind (Owocki et al. 2004). 
This uncertain mass-loss rate plagues more severely the evo- 
lution of higher-mass stars, since 15-20 M© stars stay further 
away from the Eddington limit, and overall lose little mass, 
even at solar metallicity. By what mechanism, at what rate, 
and during what phases a 100 M0 star loses mass (and an- 
gular momentum) is much less well known and this directly 
conditions the final mass and iron core mass at collapse. 

Overall, this suggests that studies of collapsar progenitors 
would benefit from a second look. Angular momentum is key 
in the current collapsar and proto-magnetar models, but there 
is a stiff requirement on the progenitor compactness to spec- 
ulate on its propensity for forming a black hole, and thus for 
producing a LGRB through one or the other channel. A ma- 
jor step forward in resolving those issues would be to con- 
duct massive-star evolution with rotation, centrifugal force, 
and magnetic fields always all the way to the formation of a 
degenerate neutronized core on the verge of collapse. This 
would allow a straightforward comparison of results between 
groups, and an easy determination of the compactness using 
GRID to test the suitability of the core for black-hole forma- 
tion. 

The ultimate check on the collapsar model requires multi- 
dimensional simulations covering the whole evolution from 
progenitor collapse, bounce, failed explosion during the 
proto-neutron star phase, formation of a black hole followed 
by the formation of a Keplerian disk, and the powering of a 
^ lOB SN ejecta. As we emphasize, black-hole formation 
is perhaps one of the most difficult steps in this sequence 
of events, and in that respect, renders the proto-magnetar 
channel quite attractive for the production of hypernovae and 
LGRBs. The diversity of LGRB/SNe, the existence of SN- 
less LGRBs and of LGRB-less hypernovae, may in fact call 
for a variety of formation channels for these rare events, in- 
cluding both collapsars and proto-magnetars. 
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